The effect of strain on the adsorption of CO on Pd(lOO) 
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We study how the binding energy, the vibrational frequencies and the adsorption isotherm of CO 
on Pd(lOO) are modified when the solid is subject to uniform strain. The parameters controlling the 
thermodynamics of adsorption (the adsorption energy, the vibrational frequencies of the adsorbed 
molecules and the interaction energy between the molecules) have been calculated by using density 
functional theory. These parameters are then used in a grand canonical Monte Carlo simulation that 
determines the CO coverage when the surface is in thermodynamic equilibrium with a CO gas, at 
a specified pressure and temperature. We find that this is substantially affected by the strain. Our 
results, along with those obtained by others, suggest that the development of "elastochemistry" , a 
study of the change in the chemical properties of a surface when subjected to strain, will lead to 
interesting and measurable results. It also suggests that differences in chemical activity between 
clusters on a support and clusters in gas phase may be partly due to the strain induced when a 
cluster is placed on the support. 



PACS: 68.10.Et, 68.10.Jy, 81.40.Jj, 82.65.My 



I. INTRODUCTION 

A number of recent experiments and calculations have 
shown that chemisorption causes surface strain and that 
straining a surface modifies the properties Cj£ [adsorbed 
molecules. Webb, Lagally and their studentallfl demon- 
strated that the structure of a sihcon surface can be 
changed by bending it. They worked with a (100) face 
of a thin silicon slab. This surface has steps separated 
by two kinds of terraces: on one, the dimer rows are 
parallel to the steps; on the other, they are perpen- 
dicular to the steps. These terraces give different low 
electron diffraction (LEED) patterns. When the slab is 
bent, the relative LEED intensities change, because the 
size of these iexraces is modified by strain. Ibach and 
his coworkersa~Ll have demonstrated a "reciprocal" phe- 
nomenon: chemisorption on the surface of a very thin 
slab will cause it to bend. If adsorption causes strain, 
then it should be possible to affect adsorption by strain- 
ing the surface. Mentzel's groupQ implanted noble gas 
ions under a metal surface, causing local strain. They 
then showed that the strained surface-has different ad- 
sorption properties. Other groupsB^Ej have strained a 
surface by growing a very thin metal film on a substrate 
of a different metal. If the two lattices are mismatched, 
but the growth is in registry, the atoms in the film are ei- 
ther stretched or compressed. This affects the chemisorp- 
tion properties of the film. In this kind of experiment the 
chemisorption properties are modified by the strain and 
also by electronic effects caused by binding to a substrate 
made of a different metal. The two effects cannot be sep- 
arated experimentally. 

Density functional calculations have been used to ex- 
plore the effect of surface strain on the prooerties of the 
adsorbates. Ratsch, Seisonen and Scheffleio have shown 
that the activation energy for diffusion of a Ag atom 
a Ag surface is affected by strain. Mattsson and Metiul 



have usad this effect to show that periodic strain on a 
surfacell3 can order nanostructures nucleated on it and 
increase their sizejjniformity. Finally, Mavrikakis, Ham- 
mer and N0rskovliZl calculated the change in the binding 
energy of O and CO and in the dissociation energy of CO 
when the Ru(OOOl) surface on which they are adsorbed 
is undeij_atrain. 

IbachtS has published a thorough review of the ef- 
fects of|-atrain in epitaxy and surface reconstruction. 
Norskovtj discusses the effect of strain in an excellent 
review of density functional studies of chemisorption sys- 
tems relevant to catalysis. 

It is not surprising that straining a surface affects the 
chemistry taking place on it. Strain changes the dis- 
tance between the surface atoms and this must change 
the properties of the adsorbates. The only question is 
whether this change can be, at least in some cases, suf- 
ficiently large to matter. If this is so, one can envision 
the development of "elastochemistry" as a new subfield 
of physical chemistry. 

In this paper we study how the binding energy, the vi- 
brational frequencies and the adsorption isotherm of CO 
on Pd(lOO) change with surface strain and coverage. We 
use generalized gradient, density functional theory for 
energy calculations and Monte-Carlo simulations and an 
analytical model for calculating the adsorption isotherm. 
The change in these quantities, caused by straining the 
surface, is sufficiently large to be measurable. In the 
case of the adsorption isotherm the coverage at a given 
gas pressure changes by almost two orders of magnitude. 
This happens because relatively small changes in adsorp- 
tion energy have large effects on the isotherm. 

To model the adsorption isotherm we need to know 
how the properties of the adsorbed molecules change with 
coverage. Ideally, we should calculate how these proper- 
ties change when we change the clustering of CO on the 
surface. One would need then the energy of all possi- 
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ble clusters. This is beyond the current capability of the 
density functional method, which uses a periodic system 
with a unit cell of limited size. For this reason we have 
adopted the following strategy. We calculate the proper- 
ties of adsorbed CO for 1/8, 1/2 and 1 monolayer. The 
one-monolayer calculation is used to determine the ad- 
sorption energy of a molecule that does not interact with 
its neighbors. The half- monolayer results provide the 
interaction energy between the next-nearest-neighbors. 
The full monolayer one is used to extract the interaction 
energy between nearest-neighbors. A Hamiltonian based 
on these energies is then used in a grand-canonical Monte 
Carlo simulation to determine the adsorption isotherm- 
We have also tested the quasi-chemical approximationES 
and found it to be inaccurate. 




FIG. 1. The square drawn with thick lines shows the su- 
percell used in low CO coverage (1/8 ML) calculations. The 
dashed lines show the supercell used for calculations in which 
CO coverage is 1 or 1/2 ML. Dark dots: Pd atoms in the top 
layer; Gray dots: Pd atoms in the subsurface layer. Black and 
open squares: the bridge sites occupied by CO. If the cover- 
age is 1/8 or 1/2 monolayer only the dark square is occupied. 
Both squares are occupied if the coverage is 1 monolayer. 



II. METHODOLOGY 

The density functional calculations repop±fld,here were 
performed with the Vienna program VASPE3'Ell This uses 
periodic boundary conditions, aj)lane-wave basis set and 
fully nonlocal Vanderbilt-typeca ultra soft psudopoten- 
tials. The generalized gradient exchange-correlation en- 
ergy is that of Wang and Perdewn^l. 

The metal is represented by a four-layer slab bounded 
by an empty region whose thickness equals that of six 
Pd layers. When the coverage p is high {p — l and 1/2), 
we use the small supercell shown by the dashed lines in 
Fig. 1; the thick lines show the supercell for the lowest 
coverage [p = 1/8). 

The CO molecules are placed on the top layer of the 
slab. The positions of the atoms in the bottom layer was 



fixed to coincide with those of the bulk Pd, calculated 
with the same program. The distance between the bot- 
tom layer and the layer next to it is also fixed to the 
bulk value. The calculated equilibrium lattice constant 
of bulk Pd is a = 3.961 A. The atoms in the other lay- 
ers are allowed to relax and take the values that give a 
minimum total energy. 

Brillouin-zone integrations have been performed on a 
grid of 9 X 9 X 1 k points for the high-coverage case and 
on a grid of 4 x 4 x 1 k points for the low-coverage case. 
The Methfessel-Paxton smearing is cr = 0.3 eV. For high 
coverage we have tested convergence by performing cal- 
culations with thicker slab and a denser k mesh. For low 
coverage, a 4 x 4 x 1 k-point mesh is the densest we can 
afford on a workstation with 1GB memory. The cut-off 
energy is 495 eV, which is a high value. 



III. THE PROPERTIES OF ADSORBED CO, AS 
A FUNCTION OF COVERAGE AND STRAIN 

A. Bulk Pd 

To determine the lattice constant of bulk Pd we used 
the supercells described in Fig. 1 and a slab having ten 
Pd layers. It is not necessary to use such a thick slab, but 
we had the results from calculations for another project. 
If the sample is unstrained, we allow the atoms to relax 
in all directions and obtain a lattice constant of 3.961 A. 
This corresponds to a distance of 1.9805 A between the 
(100) layers. We have also calculated the bulk structure 
of two samples, subject to an uniform strain of ±2% in 
a plane (the xy plane) perpendicular to the (100) direc- 
tion. In these calculations the x and y coordinates of 
the atoms are fixed, but the space between the layers is 
allowed to relax. In the sample with 2 % strain the dis- 
tance (in the z direction) between the layers is 1.9250 A; 
if the strain is —2%, this distance is 2.0457 A. Increasing 
the distance between atoms in the xy plane decreases the 
distance between the layers in the z direction. 



B. Strain and coverage dependence of adsorption 
energy 

All calculations involving adsorbed molecules are per- 
formed on a slab having 4 Pd layers, surrounded by a 
vacuum whose dimension in the z-direction is equal to 
that occupied by 6 Pd layers. The distance between the 
lowest two layers is fixed at the bulk value mentioned 
above (this distance depends on the state of strain of the 
surface). The top two layers are allowed to relax. We 
calculate the energy Es of the slabs constrained in this 
way for and ±2% strain. 

The properties of the chemisorbed CO molecules are 
calculated by placing 1 ML, 1/2 ML and 1/8 ML CO 
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on top of the Pd slab. Since these calculations are ex- 
tremely time consuming, we have investigated only the 
bridge binding sites, which have the lowest binding en- 
ergy when the surface is not under strain.cfl It is possible 
that other sites are occupied with a finite probability, 
when the system is in thermal equilibrium. Also, the CO 
molecules may prefer a different binding site, when the 
surface is under strain. These possibilities have not been 
considered in our calculations. We also ignore the fact 
that at high coverage the CO molecules might change 
their binding site. 

The energy of the CO covered slab is denoted by Etot, 
and depends on CO coverage and on surface strain. The 
adsorption energy Sad is then given by: 

£a = {Etot - Es - nEco)/n ■ 



therefore will influence the adsorption isotherm. A thor- 
ough study of these effects would have to perform a com- 
plete phonon analysis of the system. Such calculations 
are possible, in principle, but they require too much com- 
puter power, especially for the case of low coverage, to be 
attempted here. Since the adsorption isotherm is dom- 
inated by the chemisorption energy and by the lateral 
interactions, small errors in evaluating the vibrational 
energies are not likely to affect our conclusions. 



TABLE I. Ea for different coverage and strain (unit: eV) 



STRAIN 


1/8 ML 


1/2 ML 


IML 


-2 % 


-2.203 


-1.8723 


-1.20399 


% 


-2.228 


-1.918 


-1.3076 


2 % 


-2.255 


-1.944 


-1.3826 



Eco = —14.73 meV is the total energy of a free CO 
molecule, n is the number of CO molecules in a supercell 
and Eg is the energy of the Pd slab, without CO on it. 
For the supercells used in our calculations, n — 2 for 1 
ML of CO and n = 1 for 1/2 ML or 1/8 ML (see Fig. 
1). The adsorption energies calculated with this equa- 
tion, for different CO coverages and strains, are given in 
Table I. For 1/2 ML and no strain, the adsorption energy 
of —1.918 eV is consistent with that reported in ReL 04, 
which is —1.92 eV. 

The adsorption energy Ea defined above includes the 
lateral interactions between the molecules. Such a defi- 
nition is often used when the dependence of desorption 
rate on coverage is modelled. As we shall see later, this 
definition does not work well; we prefer the one given 
when we define Model II (see Section D). 



C. The dependence of the vibrational frequencies of 
CO on coverage and strain 

The dependence of the vibrational frequencies of 
chemisorbed CO on strain is of interest for two rea- 
sons. They are measurable and therefore the predictions 
made here can be tested by experiment. Furthermore, 
the change of these frequencies with strain will affect 
the chemical potential of the chemisorbed molecules and 



Because of these considerations, we calculate the vibra- 
tional frequencies by giving small displacements to the 
O and C atoms, away from their equilibrium positions, 
while keeping the metal atoms fixed. The displacements 
along the x and y directions are ±0.05A and ±0.lA, and 
those along the z direction are iO.OlA and ±0.02A. We 
calculate the forces exerted on the C and O atoms, for all 
combinations of shifted coordinates, and fit it to a linear 
form in the displacements. This allows us to calculate 
the force constants and then the vibrational frequencies. 
The results are given in Table II, for different coverages 
and strains. In Fig. 4 we show the amplitudes of four 
vibrations in the xy plane, as seen by looking down to- 
wards the surface. The highest frequency, ojg, is the CO 
stretch. In the absence of strain, at a coverage of 1/2 
ML, our result (187|8-|Cm^^) is close to that calculated by 
Eichler and HafneiO (1887 cm~^). They also report a 
frequency of 417 cm"^ for the carbon-metal stretch; our 
calculations give 341 cm^^. This discrepancy is not suf- 
ficiently large to affect the chemisorption isotherm. We 
suspect that this is due to the fact that we use more k 
points in our calculation. 

The frequency luq is easiest to measure. Our calcula- 
tions predict that either a compressive or a tensile strain, 
which changes the lattice constant by two percent, will 
increase luq by an amount measurable by infrared spec- 
troscopy. The largest change, of 29 cm~^, is at the lowest 
coverage (1/8 ML). 



TABLE II. Frequencies of CO for different coverage and strain (unit: 1/cm) 



Coverage 


STRAIN 




L02 


UJ3 


LU4 


1^5 




1/8 ML 


-2 % 


397.42 


334.23 


169.51 


49.94 


338.13 


1849.71 


1/8 ML 


% 


416.01 


346.05 


180.26 


55.86 


334.00 


1820.05 


1/8 ML 


2 % 


422.36 


327.99 


182.70 


72.12 


319.95 


1849.80 


1/2 ML 


-2 % 


407.33 


339.38 


178.37 


64.27 


395.00 


1890.61 


1/2 ML 


% 


406.60 


330.06 


170.89 


14.45 


340.63 


1878.08 


1/8 ML 


2 % 


414.24 


323.98 


182.64 


71.50 


343.95 


1883.84 


1 ML 


-2 % 


417.83 


379.86 


180.94 


56.17 


367.33 


1995.30 


1 ML 


% 


416.85 


374.93 


193.24 


77.14 


362.00 


1982.70 


1 ML 


2 % 


393.38 


356.56 


187.32 


77.00 


341.90 


1986.67 
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FIG. 2. Illustration of oscillation modes of CO in x-y di- 
rections. Black dot: C; Circle: O. 



D. The interaction between the adsorbed CO 
molecules 

To calculate the adsorption isotherm we need to know 
the energy of the interaction between molecules. A very 
precise treatment of this problem would have to deter- 
mine the energy of dimers, trimers,. . . . This would re- 
quire an enormous amount of computer power and would 
be an overkill at this stage in our research. For this rea- 
son we have studied two approximate models. 

In Model I we assume that the effect of the interac- 
tion between the adsorbates is to make the adsorption 
energy depend on coverage. This is obviously a mean 
field model; in reality the desorption energy is different, 
for different local configurations. For each strain, the ad- 
sorption energy at the three coverages used here is given 
in Table I. For other coverages we interpolate between 
these values. 

In Model II we partition the total energy into an ad- 
sorption energy and the energy of interaction between 
the adsorbed molecules. First we assume that at a cover- 
age of 1/8 ML there are no interactions between the CO 
molecules. Therefore, 

^ Et{l/%)-E,-n{ll^)Eco 
n(l/8) 

gives the adsorption energy of a CO molecule. Here 
Et{l/'&) is the total energy of the slab with 1/8 ML CO 
on it, Es is the total energy of the slab without CO, Eco 
is the energy of the gaseous CO and ri(l/8) is the number 
of CO molecules in the supercell when the CO coverage 
is 1/8 ML. 

Next, we assume that the adsorption energy of one 
molecule is independent of coverage. This means that we 
attribute the change in adsorption energy at higher cov- 
erages to the interaction between molecules. When the 
coverage is 1 /2 ML we assume that the total energy is the 
slab energy, plus the adsorption energy of the molecules 



in the supercell, plus the interaction between all the next- 
nearest-neighbors (nnn) in the supercell. Thus, we can 
calculate the interaction e„„„ between the nnn-pairs from 
the equation 

_ Et{l/2) - Es{l/2) - n{l/2)Eco - n,(l/2)g. 

Here Ua is the number of adsorbed CO molecules in 
the unit cell. The total energy of a CO monolayer can 
then be used to calculate the interaction energy between 
the nearest neighbors (nn) . We assume that the total en- 
ergy of the CO-covered surface is the energy of the slab, 
plus the adsorption energy of the CO molecules, plus the 
interaction energy of all the nn and nnn pairs in the su- 
percell. As a consequence, the interaction energy e„„ 
between the nearest neighbors is given by 

_^ Et — Eg — nnn{l)EcO — nnnn{i)£nnn — «a(l)ea 
"rm(l) 

The values of e„„ and e„„„ calculated from these equa- 
tions, for the strained and the unstrained surfaces, are 
given in Table III. 



TABLE III. Enn and Snnn dependence on strain 



Strain 






-2 % 


0.334 


0.165 


% 


0.305 


0.155 


2 % 


0.2805 


0.1555 
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We have now all the elements needed for writing an ef- 
fective Hamiltonian to be used in the statistical mechan- 
ical theory of the adsorption isotherm: the vibrational 
frequencies of the molecule, the adsorption energy, the 
energy of interaction between nearest neighbors and that 
between next- nearest neighbors. 



In equilibrium calculations we must ensure that the 
two chemical potentials use the same reference energy. 
For us, this is the electronic energy of the slab plus the 
electronic energy of a CO molecule in gas phase. The 
zero point vibrational energy of the molecule in the gas 
and on the surface is included in the chemical potential. 



IV. THE CALCULATION OF THE ADSORPTION 
ISOTHERM 

We are now ready to calculate the coverage of the ad- 
sorbed layer in equilibrium with a gas of CO molecules. 
This calculation is performed in two ways. The first 
method uses the equilibrium condition: the chemical po- 
tential of a CO molecule in gas is equal to the chemical 
potential of an adsorbed CO molecule. To calcidate the 
latter we use the quasichemical approximationcS. The 
second method, uses Grand Canonical Monte Carlo sim- 
ulations with the chemical potential of the gas calculated 
from the equations provided by Statistical Mechanics for 
ideal gases. 



A. The quasichemical approximation 

The quasichemical approximation gives for the chemi- 
cal potential of CO on Pd 



p 1 — p — a 



(1) 



where p is the coverage of CO (numbejtjDf molecules di- 
vided by the number of lattice sites) The symbol a 
stands for 



a = 2p(l~p)/(/3 + l) 



with 



/3 y^l-4p(l-p)(l-e-^""/'=«'r) . 

Since the pressure is low, we caHr-jise the ideal gas ex- 
pression for the chemical potentialE3 of the gas: 

fig = kBT\n[Pf{T)] + ^ + kgT \n{l - 6^'^°/'=-^) , (2) 



with 



ksT rnksT 2IkBT 



(3) 



In these equations, P is the pressure, m is the mass of 
CO-molecule, / is the moment of inertia of CO and loq is 
the vibration frequency of free CO, 2114 cm~^. 



By making the chemical potential for the adsorbed CO, 
[Eq. (|l|)] , equal to the chemical potential of the gas [Eq. 
(Q)], we obtain a relationship between surface coverage p 
and gas pressure P (the adsorption isotherm). Since we 
know the dependence of the parameters on strain, we can 
calculate the adsorption isotherm for a strain of ±2%, at 
300 K. 
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FIG. 3. Coverage as a function of pressure P for differ- 
ent strains. ( Solid curve: strain; dotted curve: 2% strain; 
dashed curve: -2% strain. T = 300 K. Fig. 3(a): Model I, 
Fig. 3(b): Model II.) 
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We first consider Model I, which assumes that £„„ = 
Ennn = and includcs the effect of the interactions be- 
tween the molecules in the adsorption energy Eq. The 
magnitude of Ca is given in Table I, for different coverages 
and strains. The values at other coverages are obtained 
by interpolation. For the vibrational frequencies we use 
the mean of the values given in Table II. 

The adsorption isotherm for Model I is calculated by 
setting /ig given by Eq. (||) equal to /la given by Eq. (|^), 
with e„„ = and ea replaced with e^. This gives us 
an equation connecting the gas pressure to the coverage 
(the adsorption isotherm). The dependence of coverage 
on the gas pressure obtained with this method is plotted 
in Fig. 3(a), for and ±2 % strains, at T = 300 K. 

In the next calculation we use Model II, which assumes 
that the adsorption energy ea is independent of coverage, 
that the nearest -neighbors interact with the energy e^n^ 
and that e„„„ = 0. The isotherm is obtained by making 
^^a = jJ'g- The resulting adsorption isotherm is plotted in 
Fig. 3(b). 

It is clear that in both models the strain has a very 
large effect on chemisorption. For example, the Fig. 3(b) 
shows that at 10~^ Torr the CO coverage on the un- 
strained surface is ~ 0.64 ML. Stretching the metal to 
increase the lattice constant by 2% changes the coverage 
to ^ 0.92 ML. Another way to look at this is to note that 
a coverage of 0.7 is achieved at roughly ~ 1.2xl0~® Torr 
if the lattice constant is increased by 2%, at ~ 1.2 x 10^^ 
Torr if the surface is not strained, and at ^ 10"^ Torr 
if the lattice constant is decreased by 2%. These effects 
are large enough to be easily detected experimentally, if 
one could find a convenient way to strain the surface in 
ultra-high vacuum. 



B. Monte Carlo simulation 



for 2 % strain. The effect is much stronger if the quasi- 
chemical approximation is used (see Fig. 3 (b)), but we 
attribute this to inaccuracies in Model II. The discrep- 
ancy between this model and the Monte Carlo simulation 
is so great that we have to discount Model II as insuf- 
ficiently accurate. We suspect that this is due to the 
neglect of the repulsive, next-nearest-neighbor interac- 
tions, which makes the coverage given by Model II larger 
than in reality (here the Monte Carlo simulation is the 
"reality" ) . Model I is so far off from both the quasichem- 
ical approximation and the Monte Carlo simulation that 
it should be completely discarded. This is a pity since, 
had it been correct, this model would have allowed a very 
simple interpretation of the experiments. 



Monte Carlo Simulation 
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FIG. 4. MC simulation of coverage as a function of pres- 
sure P for different strains under the assumption of no nn 
interaction between COs. T = 300 K. 



We use grand-canonical Monte Carlo (MC) calcula- 
tions to determine the coverage of CO on the Pd(lOO) 
surface, at given gas pressure and temperature. The 
methodology is presented well in an excellent book by 
Frenkel and SmitE3 so there is no need to discuss it here. 
The Hamiltonian consists of the adsorption energy 
multiplied by the number of adsorbed molecules, plus the 
pairwise interaction energy e„„ multiplied by the number 
distinct nearest-neigbor pairs, plus the energy e„„„ mul- 
tiplied by the number of distinct next-nearest-neighbor 
pairs, plus the adsorption energy, plus the vibrational 
energies of the adsorbate. The chemical potential of the 
gas is calculated from Eq. (||). We use 1000 x 1000 lattice 
sites. 

The results of the simulations are plotted in Fig. 4 for 
T ~ 300 K. As in the calculations performed in the pre- 
vious section, we find that strain substantially affects the 
surface coverage. For example, at a pressure of 10~^ Torr 
and a temperature of 300 K, the coverage is ~ 0.70 for 
—2 % strain, ~ 0.66 in the absence of strain, and ~ 0.57 



V. SUMMARY 

We have studied how uniform strain in the plane paral- 
lel to the (100) surface of palladium affects the properties 
of chemisorbed CO. The effect on the vibrational frequen- 
cies is small: it is comparable to or less than that ob- 
served for the molecule located in different binding sites. 
The shift in the stretch frequency of CO is sufficiently 
large to be detected by infrared spectroscopy. The effect 
on the binding energy at low coverage is larger. We found 
the weakest binding when the surface is compressed (i.e. 
for —2 % strain). Releasing the compressive strain raises 
the binding energy by 0.025 eV. Going to a tensile strain 
of 2% increases the binding energy by another 0.027 eV. 
From the equilibrium condition ( fig given by Eq.(^ is 
equal to fiad given by Eq. (|l|) ) we see that the logarithm of 
the equilibrium pressure is proportional to the adsorption 
energy. This means that small changes in the adsorption 
energy lead to large changes in the equilibrium pressure. 
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While the quasi-chemical approximation is not quantita- 
tively accurate, this particular prediction is reliable. It 
is indeed confirmed by the Monte Carlo simulations. 

It is interesting to speculate on the broadej- conse- 
quences of these findings. Ibach's experimcntscl^El showed 
that chcmisorption induces surface strain, which results 
in the bending of the thin slab on which the adsorption 
was performed. This means that in all equilibrium stud- 
ies we ought tp_take strain as an additional thermody- 
namic variableEZi and stress as its conjugate variable. It 
follows that all equilibria will depend on the strain state 
of the surface. Moreover, the outcome of transforma- 
tions performed at constant strain will be different from 
that of transformations taking place at constant stress, 
just as transformations at constant volume differ from 
those at constant pressure. This will affect the adsorption 
isotherm and chemical equilibrium at the surface. Since 
strain affects equilibrium, detailed balance tells us that 
it must also affect kinetics. The only question is whether 
the modifications caused by strain are large enough to 
make it necessary to include them in a thermodynamic 
analysis of surface phenomena. The present work sug- 
gests that they are, and therefore they deserve further 
study. 

Besides using simple mechanical means, one could in- 
duce strain by passing sound through the solid. In par- 
ticular, the excitation of the Rayleigh mode would seem 
most effective, since it affects the surface more than 
other modes. Based on our calculations one would ex- 
pect desorption to occur from the regions compressed by 
the sound wave. Moreover, chemical equilibrium on the 
surface is-iikely to be affected. Indeed, Krishner and 
LichtmanE3 have shown that passing sound through a 
surface causes desptption of the molecules adsorbed on 
it. Inoue et aZ.t5L3 have demonstrated that sound af- 
fects the yield of catalvtic-reactions. Prompted by these 
experiments King et a/oO have performed careful ultra- 
high vacuum experiments in which the interplay between 
sound and chemisorbed molecules can be studied in a 
clean and controlled environment. They found that ex- 
posing a Pt surface to sound can change the rate of CO 
oxidation. By using photo-electron rpirroscopy, Kelling, 
Cerasari, Rotermundt, Ertl and King^ have shown that 
the reaction is affected mostly by controlling CO adsorp- 
tion at given CO pressure. 

In most "sonochemistijy" , sound acts by heating and 
compressing the systemcH. This is not the case in the sur- 
face science experiments jpeationed above, which have 
ruled out heating effcctsEjEa. The effect of sound on 
chemical reaction is very hard to explain by conven- 
tional arguments: sound frequency is very low and ex- 
cites phonons of long wavelength with very low energy. 
It is unlikely that these can affect the rate of chemical 
processes. For this reason we speculate that perhaps it 
is the strain induced by the sound wave that plays a role 
in the process. The fact that the presence of sound does 
not affect the activation energy of the oxidation reaction 
but influences the CO coverage supports t his view. 



We expect that strain effects also play a role in sup- 
ported catalysts, especially those using small clusters. 
When such a cluster is deposited on a support it suffers 
two modifications: a charge rearrangement and a distor- 
tion of its geometry (strain) . It is therefore possible that 
a part of the change in the catalytic activity cause by 
depositing clusters on a inert support comes from the 
strain. 
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